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Abstract 

We formulate the Hubbard model for the simple cubic lattice in the 
representation of interacting dimers applying the exact solution of the 
dimer problem. By eliminating from the considerations unoccupied 
dimer energy levels in the large U limit (it is the only assumption) we 
analytically derive the Hubbard Hamiltonian for the dimer (analogous 
to the well-known t — J model), as well as, the Hubbard Hamiltonian 
for the crystal as a whole by means of the projection technique. Using 
this approach we can better visualize the complexity of the model, 
so deeply hidden in its original form. The resulting Hamiltonian is a 
mixture of many multiple ferromagnetic, antiferromagnetic and more 
exotic interactions competing one with another. The interplay be- 
tween different competitive interactions has a decisive influence on 
the resulting thermodynamic properties of the model, depending on 
temperature, model parameters and assumed average number of elec- 
trons per lattice site. A simplified form of the derived Hamiltonian 
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can be obtained using additionally Taylor expansion with respect to 
x = jj (i-hopping integral between nearest neighbours, ?7-Coulomb 
repulsion). As an example, we present the expansion including all 

,2 

terms proportional to t and to jj and we reproduce the exact form of 
the Hubbard Hamiltonian in the limit U — > oo. 

The nonperturbative approach, presented in this paper, can, in 
principle, be applied to clusters of any size, as well as, to another 
types of model Hamiltonians. 

1 Introduction 

The single-band Hubbard model, Ref. [1], plays in the solid state physics 
a similar principal role as the hydrogen atom in the atomic physics. This 
explains a continuous interest in its properties. The Hubbard Hamiltonian 

reads _ 

i,j,a i 

Here c ija (cf a ) are annihilation (creation) operators of an electron with spin 
a =f, J, in the Wannier representation at the lattice site Rj and n ita = cf a Ci ta . 
Moreover, tij is the hopping integral between different lattice sites i and j 
(U,i = 0) and U is the intrasite Coulomb repulsion. The Bloch conduction 
band energy is given by 

,< »„' ' k H H • (2) 

i-j 

In the following we restrict ourselves to the simple cubic (sc) lattice and 
assume that 

f — t i,j-nearest neighbours 
% ' j \ otherwise. 

Then 

£k = — 2t(cos k x a + cos k y a + cos k z a). (3) 

The hopping parameter t is simply related to the bandwidth W of the 
Bloch band (2), e. g. W = 12t for the sc lattice. The interplay between the 
two model parameters, W and U, is decisive for the properties of the model 
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resulting in strong electron correlations, leading to band magnetism (see e.g. 
Refs [2], [3] for a review), insulator-to- metal transition (see e.g. Refs [3], [4] 
and papers cited therein) and high-T c superconductivity (negative [/-model, 
see e.g. Ref. [5]). The Hubbard model very often plays also a role of a 
submodel for many other more complicated models (as e.g. Anderson model 
(see Ref. [6]), s-f model (see Refs [7,8]) and so on). Especially interesting, but 
difficult to handle are the properties of the Hubbard model in such a regime of 
the model parameters where the bandwidth W is comparable to the Coulomb 
repulsion U. In a large number of papers [9-30] the authors tried to solve this 
model using many sophisticated methods. The exact solution, however, does 
not exist till now. Many authors tried to change the situation in this field 
by introducing the expansion parameter x = jj (i C 1). This idea (cf Refs 
[31,32]) consists in replacing "difficult physics" connected with the model by 
"difficult mathematics" obtained by a laborious expansion with respect to x. 
Different methods connected with this problem have been applied as e.g. the 
perturbation expansion (see Refs [18], [21]), canonical transformation (see 
Refs [9], [12], [17], [22], [23]) or ab initio derivations (see Refs [25], [28-30]). 
Most of the methods, leading to the t — J model (or generalized t — J model) 
are also summarized in Refs [27], [4]. 

The goal of the present paper is just to show that we can take another, 
nonperturbative way. In the first step we divide the crystal lattice into a set of 
interacting dimers. In other words, we can rewrite the Hubbard Hamiltonian 
(1) for the sc lattice (see Fig. 1) in the equivalent form Q 

I h° 

where 

H i = -tJ2( c i,i,* Cl > 2 >° + c t,2,a c i,i,°) + U(ni M ni M + n Ji2)T n Ji24 ). (5) 

CF 

-'-The dimer Fourier transformation 

k k 

where N is the number of lattice points, applied to (4), gives the well known result 

H = J2 £knu,a + f E C k+q,T C k-T C k'-q,l Ck '.i Wlth £k S iV<3n ^ ( 3 )' 
k,(T k,k',q 
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The indices / and J enumerate the dimers and Hf is the dimer Hamiltonian. 
The second term in (4) describes the hopping between nearest dimers in 
the ^-direction, the third one represents the hopping between nearest dimers 
((?/, z)-plane) and between different dimer planes (see Fig. f). 




Figure 1: A plane of the Hubbard dimers. The crystal is considered as a collection of 
parallel lying dimer planes. 

The equivalent form of the Hubbard Hamiltonian (4) is especially suit- 
able because we can apply in the following the exact solution of the dimer 
problem (5). In the next step we express the construction operators a a (a+) 
as linear combinations of the transition operators between different dimer 
states. This, in turn, allows to find the exact dimer representation of the 
whole Hamiltonian (4). The space of the dimer eigenvectors consists, how- 
ever, of two subspaces. One of them corresponds to the lowest lying energy 
levels, the second one contains the levels with energies which in the large 
U limit take on large, positive values. These levels in a reasonable temper- 
ature range cannot be occupied by electrons and therefore we can exclude 
them from further considerations. It is interesting to note that this approch, 
without any additional assumptions, applied to the dimer Hamiltonian (5) 
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produces the analogy of the well known t — J model (cf Refs [31], [32]). A 
similar approach can be applied to the Hamiltonian (4), describing the crys- 
tal as a whole. By eliminating from the considerations the unoccupied dimer 
energy levels in the large U limit (it is the only assumption) with the use 
of the projection technique we can find in a straightforward way the final 
form of the Hamiltonian in this limit without using perturbation expansion 
or canonical transformation. The resulting Hamiltonian, obtained in this 
way, is very complicated. It, however, explicitly shows all possible magnetic, 
nonmagnetic or more complex competitive interaction processes very deeply 
hidden in the original Hamiltonian written in the site representation (1). It 
is the aim of this paper just to reveal these important but normally invisible 
elementary interactions. One important advantage might be that a given 
approach to the unsolvable Hubbard problem can be tested with respect to 
the types of neglected interaction processes. Besides, the new and straight- 
forward method, presented in this paper, can easily be adopted to clusters 
of any size and also to another types of model Hamiltonians. 

The paper is organized as follows. In Sec. 2 we find the exact solution 
of the Hubbard dimer (5) and give the exact expressions for the annihila- 
tion operators cj j i(2), ( j in the dimer representation. In Sec. 3 we derive the 
Hubbard dimer Hamilonian (5) in the large U limit. With the use of the 
projection technique onto the lowest lying dimer states we derive in Sec. 4 
the Hubbard Hamiltonian for the crystal in this limit (central formula of this 
paper). A simplified version of the derived Hamiltonian with the use of the 
Taylor expansion with respect to x = ^ (i C 1) and the case U — > oo is 
discussed in Sec. 5. 

2 Exact solution of the Hubbard dimer 

The eigenvalues and eigenvectors of the dimer Hamiltonian (5) can be found 
by using a standard procedure (cf Refs [33-35]). We start with the vectors 
|rii,T, ni,4.; ri2,T, fi2,4.) ip>i,c — 0, 1; i — 1, 2; a =1,1) which constitute the Fock 
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basis of the single-dimer space of states: 



o) = 


= |0,0;0,0), 


|21> = 


|1,1;0,0), 


|31) 


= |0,1;1,1>, 






|22) = 


|1,0;1,0), 


|32) 


= 1,0; 1, 1), 


11) 


= |1,0;0,0>, 


23) = 


|1,0;0,1), 


|33) 


= lM;0,i), 


12) 


= |0,1;0,0>, 


24) = 


|0,1;1,0), 


|34) 


= |i,i;i,0>, 


13) 


= |0,0;1,0), 


25) = 


|0,1;0,1), 






14) 


= |0,0;0,1), 


26) = 


|0,0;1,1), 


|4) = 


= |i,i;i,i>- 



(6) 



Starting with the vectors (6) we easily get the eigenvalues E a and the eigen- 
vectors \E a ) of the Hubbard dimer (5). Here, we only mention that the 
space of 16 eigenvectors \E a ) can be devided into some subspaces num- 
bered by n = J2 n i,v The subspace, belonging to n = and n = 4 is 

i,cr 

1-dimensional (E = 0, \E ) = |0); E A = 2U, \E A ) = |4)). There are, how- 
ever, two 2-dimensional subspaces corresponding to n — 1 (as e.g. En = —t, 
\En) = -^=(|11) + |13)), etc.) and two 2-dimensional subspaces correspond- 
ing to n = 3 (as e.g. E 31 = t + U, \E 31 ) = ^(|31) + |33)), etc.). The 
subspace, belonging to n = 2 consists of one 4-dimensional subspace (as 
e.g. E 2 i = 0, I-E21) = 75(123) + 1 24) ) , etc.) and two 1-dimensional sub- 
spaces (£"25 = E 2 q = 0, I-E25) = 1 22) , 1^26) — 1 25) ) . The complete set of the 
eigenvalues E a and eigenvectors \E a ) of the Hubbard dimer is given in the 
Appendix A. It allows to express the dimer Hamiltonian (5) in the equivalent 
form 

H d = J2E a \E a )(E a \. (7) 

a 

The next important step in our calculations is the possibility to express 
the annihilation operators cu2\ a as linear combinations of transition oper- 
ators between dimer states P a ^ = \E a )(Ep\. This procedure can easily be 
performed when acting with the annihilation operators on the basis vectors 
(6) and using the reciprocal relations to (A.l). In this way we can obtain 
the dimer representation of the annihilation (creation) operators, given in 
Appedix B. After this operation we can insert so prepared c/ j i(2), - ( c / u 2 ) a) 
into the Hubbard Hamiltonian (4) to obtain the Hubbard model in the dimer 
representation for the sc lattice. This representation will be used later to de- 
rive the Hubbard Hamiltonian in the large U limit (see Sec. 4). 
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3 Hubbard dimer for large U 



Looking at the eigenvalues (A.l) of the Hubbard dimer it is easy to see that 
for large U (U > t) the energies E a — E22, E23, -E31, -E32, -£"33, -£"34 and E^ 
take on large, positive values, producing in the partition function the terms 
which can practically be neglected. It means that the mentioned energies are 
not occupied in the reasonable temperature range (1 eV ~ 11604.5 K) and 
can be excluded from our considerations. Therefore the dimer Hamiltonian, 
given by (7), reduces to 



H d = -tP 11 , 11 + tP 12 ,12 - tP 13 , 13 + tP UM + \ -C + -J P 24 ,24. (8) 

To bring this expression into a compact (second quanization) form we intro- 
duce the Hubbard operators 



and spin operators 



(Li,cr Q,<t(1 ^i,—a)i (9) 

bi, a = Ci^Tli-a (10) 



si = ^Kt-^) = ^k t -<j, (ii) 

s t = c ti c i,i = a ti a i,i> ( 12 ) 

^ = c+Q,T = <|^,T (13) 

where n a ia = a+ a a ij<T (i = 1,2; a =|, |). 

With the use of (9)- (13) the Hamiltonian of the Hubbard dimer (8) for 
large U can be presented in the form (we introduce the omitted earlier dimer 
index /) 
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a 



H / = • & 1,2 7 J 

'4i\2 4 



^ 2 (i- V 1 + (#) 2 ) 



+ b t,2,i a U,i a i,i,l h iA,l) ( 14 ) 
+ (1 - <,i)<2 + (1 - n?>J,i - <i< 2 ] 



J2 Yl ["U* + ^a.^wl 



^ 2 a a=l,2 



where raj£ = n l«,^ = K^i^a = n hoi ^n hoi ^ c (a = 1,2), 

a — 1 when a = 2 and a = 2 when a = 1. It is very important to stress that 
the formula (14) has been obtained in a nonperturbative way, starting from 
the exact form of the dimer Hamiltonian (see (5) or (7)), and excluding from 
the considerations unoccupied dimer energy levels in the large U limit. Let 
us note that using the Taylor expansion in (14) with respect to x — jj and 

retaining the terms proportional to jj we obtain 



Hi = -tEKi^A. + <2,^/,i,.] + -plS i,i ■ Sj, 2 - (15) 

cr 

The formula (15) is the well-known t — J model for the Hubbard dimer where 
the first part in (15), similarly to (14), represents the exact form of the dimer 
Hamiltonian (5) in the limit U — > oo. 

The same result (14) can also be obtained applying a more general ap- 
proach which will be used later to derive the Hubbard Hamiltonian for large 
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U in the case of a crystal. Let us note that after elimination of the unoc- 
cupied levels the subspace of the eigenvectors of the Hubbard dimer (A.l) 
for large U consists of the following eigenvectors: l-Eo), l-Eii); \Ei2), I-E13), 
\Eu), I-E21), I-E24), I-E25), and l-E^)- It means that we can define a projection 
operator onto this subspace 



(16) 



P _ p(-o , P (n , P v) , P (-o , p(^) 

— -"0,0 ^11, 11 ~r ^12, 12 ~r ^13,13 "r ^14,14 

, P (n , p(-f) , pco , P oo 

"l - 21,21 "r ^24,24 "r -"25,25 "r -"26,26 

which, in the second quantization form, reads 

p I= l-K^+fiy + M,^ 

+ 2 (^i,T<ia a J,2,^/,2,T + ^ 2 , T <2,| a /,i,A,i,t) (17) 
+ (1 - n?,iK, 2 + (1 - n? )2 Ri - <i< 2 ] 

+ rr /TT71^2 S ^ ( a i>,cA ',a,<r + H,a,cr a I ,a,a) ■ 
u \ ,1 +(-U> a a=l,2 

Now, the Hamiltonian (14) can also be obtained from (5) (or (7)) with the 
use of the projection operator (16) or (17) using the relation 

Hi = PjHfPj (18) 
and applying a straightforward but laborious algebraic calculation. 

4 Hubbard model for large U 

The Hamiltonian of the whole crystal (4) can be expressed (similar to (7)) 
in the form 

H = ^2e i \E 1 )(E 1 \ (19) 

7 
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with unknown energies E 1 and eigenvectors \E 1 ). We can, however, expand 
the eigenvectors |i? 7 ) in the series of the dimer eigenvectors (see (A.l)) 

l^> = E c l,,-y M \ E ^)-\E lM ) (20) 

71, -,7m 

assuming that the crystal consists of M dimers. Using (19) and (20) we 
obtain 

" = E^ E E C 7i, -,7m C 7i , -,7m I ^71 ) • • I E lM ) I " ( E 1' M I • ( 21 ) 
7 71j--)7m7i,--,7m 

It is clear that to obtain the Hubbard Hamiltonian for large U we have 
to project (21) onto the subspace of the lowest lying dimer states with the 
use of the projection operator 

P = P 1 P 2 ...P M (22) 

where Pj is given by (16) or (17). In analogy to (14) we denote the Hubbard 
Hamiltonian in the large U limit by H. Similar to (18) we write 

H = PHP (23) 

and instead of the form (21) for the Hamiltonian H we can use (4). Taking 
into account that P 2 = P (Pj = P u \P U Pj\ = 0) we obtain 



H = P[J2 H r -tJ2 «2,<x C m,l,- + 4+l,l,a C I,2,a) 

1 _ t'l = (24) 

~t E (4,l,a C JM + 4,2,a C J,2,a)] = PH 

where is given by (14) and (a — 1, 2) 

Cl,a,a — Cl,a,<jPl, (25) 
4,a,« = Picf^r- (26) 

Applying the projection operator (16) or (17) to (B.l) - (B.4) and intro- 
ducing Hubbard- and spin operators (9) - (13) we obtain 

c/,i )T = + P [5| i2 aj,i )T + Sj i2 a Ixl ] - 5 [Sj^^ + Sj^ai^] , (27) 
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c/,u = ~ P [S", 2 ai,i,i ~ ^ 2 a/,i )T ] + 5 [5f ;1 a 7 , 2a - S^a/,^] (28) 

where the corresponding expressions for ci,i,a ( a — 1\ I) can easily be 
obtained by changing the internal dimer index 1 <^ 2 in (27) and (28). 

The new operators a (a = 1,2; cr =T, I) are introduced to obtain a 
relatively compact form of (27) and (28). They are defined as follows 



a 



where 



[ , nf ] (29) 

+<-> 0/,a, CT + a l t a-a a I,a-aOl,a,a + aj^-^- 



<W = a w ( 1 - f - ^ ) , (30) 
&W = W (31) 



and 




/5 = 7( 1 -^=^] ) (32) 

(33) 

To write down the explicit form of the Hubbard Hamiltonian H (see (24)) 
we have to insert (27) and (28) into H. This operation leads, however, to 
a very complicated form of H, given in the Appendix C (a simplified form 
of this Hamiltonian is discussed in the next Section). Here again the only 
assumption made to derive H was the reduction of the whole dimer space 
(A.l) to the subspace of the dimer eigenvectors \Eq), I-E12), \Eiz), 

\E 14 ), \E 2 i), I-E24), I-E25) and \E 2 g), corresponding to the lowest lying dimer 
energy levels because in the large U limit only these levels can be occupied. 
The Hamiltonian H, obtained in this way, contains many competing, mag- 
netic, nonmagnetic and more complex interactions. Among them we can 
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find a direct antiferromagnetic interaction generated by the term Si,i ■ Si,2 
(Heisenberg exchange interaction) multiplied by the positive coupling con- 
stant. Such a term appears in H 1 (see (C.6) and (14)). Inside of H (see 
(C.6) and (C.2) - (C.5)) a kind of ferromagnetic interactions between spins 
from different dimers, represented by (C.3), appears within the terms pro- 
portional to (3 2 and 5 2 (negative coupling constants). The antiferromagnetic 
interactions, however, appear again in terms proportional to (35. There are 
also many other magnetic, more exotic interactions, represented by (C.2), 
(C.4) and (C.5), entering into (C.6). The situation is, however, much more 
complicated when we consider the total Hamiltonian H (see (24)) in the large 

U limit. H differs from H by the multiplicative factor P (a product of the 
projection operators Pj (see (22) and (17)), standing on the left. Inside of 
each Pj the mentioned antiferromagnetic interaction also appears (see (17)). 
In other words, the total Hamiltonian H, we are interested in, is actually a 
sum of the products of many competitive ferromagnetic, antiferromagnetic 
and more exotic interactions. The thermodynamic properties of the system, 
described by the Hamiltoniam H (24), are then a result of the competition 
between all of them. Which interaction wins in such a competition it cer- 
tainly depends on temperature, model parameters (t, U) and on the average 
number of electrons per lattice site, determining the chemical potential of 
the system. 

The formalism presented in this paper is also applicable to a more compli- 
cated decomposition of the Hubbard Hamiltonian (1) into a set od interacting 
clusters consisting e.g. of one central atom and z its nearest neighbours. We, 
however, know (see e.g. Refs [36]- [39] and papers cited therein) that, un- 
fortunately, the mathematical problems in this case exponentially grows up 
with the size of the cluster. 

5 Taylor expansion 

The complicated form of the Hubbard Hamiltonian H in the large U limit 
(see (24), (C.6)) where P is given by (22) (see also (17)) can essentially be 
reduced when applying the Taylor expansion with respect to the parameter 
i = ^ < 1. To do it we have to expand all the coefficients in (14), (17) 
and (C.6) including also f3 and 5 (see (32), (33)). Such an expansion can 
be performed to any power of x, however, the most simple form we obtain 
when we restrict ourselves to the linear approximation, resulting in the terms 
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proportional to t and jj. The accuracy of this expansion can easily be verified 
assumming e.g. a typical value of the ratio = \ ( or less). Because the 
bandwidth of the conduction band for the sc lattice is W = 12t, it results 
in a small value of the expansion parameter x = jj = ^ in this case. It. 
however, means that the linear approximation is quite reasonable because all 
higher terms in the expansion, proportional to x n (n = 2,3, ...), produce 60 
times smaller contribution. 

To present the results of the Taylor expansion including all the terms 
proportional to t and jj let us first define several auxiliary quantities 

M 



pW = i^F«.....p<l\pf)FW.....p« (35) 

1=1 



1=1 

M 



where 



^^^Ki+^y+^A ( 36 ) 

2 

P i 2) = E E (°W 6 W + bta^i,^) (37) 

a a=l 



and (a = 1, 2) 



"J.a.o = "/ a.a [ 1 ~ ) ■ ( 3S ) 



Ul,a,a = kl,a,c + a l,a-a a I,*,-abl, a ,a + O/^ct^- (39) 

where 6/ )Q , )(7 is given by (31). 

The Hamiltonian if (C.6), including all the terms proportional to t and 
jj : takes on the form 

T = + 17 (2) (40) 
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where 



and 



I, a 



- l Yl [ a h,<T a W,° + a i+l,l,a a I,2,a] (41) 



I,a 

2 



I^J,cr a=l 



+ a /+l,l,<7 a J,2,ff + a /+l,l, CT a -T,2,o-J 

~77 E Z^[ a W a ^ + a w aj,a, CT ] (42) 

I^J,a a=l 
2t 2 — > 

7 



+ 5 7+1,1 • (~| J+l,2;/,2 + "?/,2;7+l,2)] 



2t 2 — > 



+ ^7, 2 -(^W,2 + ^ > J,2;/,l)]- 



The operators s z,± and //, ,„ in (42) retain their forms introduced in 

(C.l) but a j ^ o .(a / f ' a ) in (CI) should actually be replaced by ar^aifif ^ a ), 

defined by (38). The total Hamiltonian H in the large U limit (24) including 
all terms proportional to t and jj can thus be written in the form 

h = pm(S {1) + 1 (2) ) + p^t (1) = p^m {1) + (p«i7 (2) + p^m {1) ). (43) 
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x=(l) 

The first part, P^'H , contains the terms proportional to t whereas the 
term in the parentheses is proportional to jj (see (35) and (42)). In the 

=(2) 

expression for H (see (42)) the first term describes the antiferromagnetic, 
Heisenberg intradimer interaction (see also the second term in (15)). The 
apearence of this term may suggest that such an interaction should also arise 
between different dimers. This is of course the case. However, because of 
applied procedure such terms do not appear explicitly. In the last analysis 
our method treats all interactions within and between dimers, on the same 
quality level, i.e. all interactions are taken into account. The magnetic 
interdimer interactions, represented by the fourth and fifth term in (42), have 
formally the same structure as the Heisenberg interactions but instead of the 
scalar products of spin operators there are the products of spin operators 
and "hopping spin" operators (defined in C.l). All the terms presented in 

(42) are correct because they originate from the exact decomposition of the 
Hubbard Hamiltonian (1) into a set of interacting dimers (4) where each 
dimer problem has been exactly solved (exact dimer representation of the 
construction operators after applying the projection procedure, given by (27) 
and (28)). The total Hamiltonian H (43), we are interested in, is much more 

=(1) =(2) 

complicated than H and H alone (see (41)- (43)) because of the presence 
of the projection operators and P^ (see (34)- (37)) in the expression 

(43) . 

It is also interesting to see what happens in the special case when taking 
the limit U — > oo. The second term in the parentheses of Eq.(43) vanishes 
in this limit (cf (35) and (42)). Besides, each lattice site cannot be at the 
same time occupied by two electrons what is equivalent to the assumption 
(a =1,2) 



n. 



and (cf (38)) 



b I,a,a b I,a,a = ™/,a )( T™/,a -a = 0, (44) 

CT 

ai jCe , a = a/ iQ)CT . (46) 
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The total Hamiltonian (43) is thus given by a simple formula 
H = P^T = -t £ [a+^ar^ + aj^a^] 



I, a 



~t J2 [ a t,2,a a I+l,l,v + a t+l,l,a a I,2,o] (47) 
2 

~t a t,a,a a J,a,v 
IjLJ^a a=l 

Going back to the original lattice (cf (4), (5) and (1)) it is easy to see that 
the exact Hubbard Hamiltonian in the limit U — > oo takes on the form 

H=-tJ24,*H° (48) 

where i and j (as before) number the lattice points and Oj i(7 (a^ a ) are the 
Hubbard operators, defined by (9). 



6 Conclusions 

Using a new, nonperturbative approach, basing on the equivalent form of the 
Hubbard Hamiltonian, represented by the collection of interacting dimers (4) 
where each dimer problem has been exactly solved, we have expressed the 
annihilation (creation) operators (B.1)-(B.4) as linear combinations of tran- 
sition operators between different dimer states (dimer representation). This 
method made it possible to exclude from the considerations the unoccupied 
dimer energy levels in the large U limit by means of the projection technique 
resulting in the final Hamiltonian for the dimer itself (14), as well as, for the 
crystal as a whole (see (23), (24) and (C.6)). It is important to stress that the 
elimination of the unoccupied dimer energy levels was the only assumption 
to derive the Hubbard Hamiltonian (24) in the large U limit. Therefore we 
can be sure that expanding (24) with respect to x = jj (Taylor expansion) 
we obtain absolutely all terms proportional to x n (n = 1,2,...). In other 
words all the coefficients proportional to x n are easy to control what is not 
always the case when using another methods. The final form of the obtained 
Hamiltonian in the large U limit (see (23), (24) and (C.6)) visualizes high 
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degree of complexity of the model, deeply hidden in its original version, form- 
ing a mixture of many multiple ferromagnetic, antiferromagnetic and more 
complex interactions competing one with another. This fact seems to be 
decisive for our final conclusion. Because we are still dealing with the ap- 
proximate solutions of the model (the exact solution does not exist till now) it 
may happen that we underestimate in this way some important interactions 
and overestimate the others. It is the reason why the resulting thermody- 
namic properties of the Hubbard model, obtained in an approximate way, so 
strongly depend on the quality of applied approximations. 

We note in passing that the exact dimer solution may serve as a novel 
alloy analogy for the Hubbard model, which could be treated by coherent 
potential approximation (CPA). It is well-known that the standard alloy 
analogy, based on the atomic limit, does not allow for ferromagnetism in 
the Hubbard model. This may change by application of the dimer solution 
which already accounts for a restricted hopping of the band electrons. A 
corresponding study is in preparation. 

Another example of a dimer approach to Hubbard-like models is the bond 
operator theory as an extension of the slave bosonic and fermionic operators 
(see Refs [40], [41] and papers cited therein). 
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Appendix A 



The exact solution of the dimer eigenvalue problem (5) reads: 

E = 0; \Eq) = \0), 



E n = -t; |S 11 > = ^(|11> + |13», 

En = t; \E 12 ) 

E13 = —t; \E 13 ) 

E u = t] \Ei4) 



V2 

1 



(111) 



|13», 

^(|12> + |14», 
^(|12>-|14», 



E21 
E22 

E23 
E24 
E25 
E26 



where 



0; 
U; 

C+f, 

-c + f; 

0; 

0; 



E31 
E32 
E33 
E34 



E21) = ^(|23) 
E22) = ^(|21) 
E23) = ai(|21) - 
E 2i ) = a 2 (|21) - 
£25) = |22), 
E 2 6) = |25), 



■|24», 
■|26», 
|26» - 
|26» ■+ 



a 2 (|23) 
«i(|23) 



|24», 
124)), 



(A.l) 



t + U; 
-t + U; 
t + U; 
-t + U: 



\E 31 ) 
\E 32 ) 
\Eaa) 

\E3A) 



^(|31) + |33», 
^(|31>-|33», 
^(|32> + |34», 
^(|32>-|34», 



E A = 2U- \E A ) = \A) 




a2 = \v-^c- 



(A.2) 
(A.3) 
(A.4) 
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Appendix B 



The exact dimer representation of the construction operators is given by 
the following expressions: 



C l,t = 7f (-fo,ll + ^0,12) + 75^11,25 - 75 -Pl2,25 (B.l) 
+ | (-Pl3,21 + ^13,22) + 75 (^-Pl3,23 + ^13^4) 
— g (-Pl4,21 — Pl4,22) + 75 (aPl4,23 ~ ^14,24) 
+ | (-^21,33 — -P21,34) — 2 (-^22,33 + -^22,34) 
+ 75 ( a -P23,33 + ^23,34) — 75 (^24,33 — aP 2 4,34) 
+ 72 ("^26,31 - ^26,32) + ^31,4 + 75^32,4, 



Cl,i = 75 (-^0,13 + ^0,14) + 75-Pl3,26 _ 72=-Pl4,26 (B-2) 
+ | (-Pll,21 — -Pll,22) ~~ 75 (&-Pll,23 + O-Pllfli) 
— \ (-Pl2,21 + -Pl2,22) — 7J (0-Pl2,23 — ^-Pl2,24) 
~| (-^21,31 — ^21,32) — I ("^22,31 + -^22,32) 
+ 75 (0-f23,31 + ^23,32) — 75 (^24,31 — ^24,32) 
-75 (-^25,33 - ^25,34) - 71-^33,4 ~ 75^34,4, 



c 2,t - 75 (-Po.ii _ ^0,12) - 71 -Pli,25 - 75^12,25 (B.3) 

— I (-Pl3,21 + -Pl3,22) + 75 (^-Pl3,23 + 0-^13,24) 

— 2 (-Pl4,21 — ^14,22) — 75 ( a -Pl4,23 ~ bPu^A) 

— \ (^21,33 + -^21,34) + \ (-^22,33 _ -^22,34) 

+ 75 ( 0-^23,33 — ^23,34) — 75 (^24,33 + oP 2 4,34) 
-75 (^26,31 + ^26,32) + 75-^31,4 - 75-^32,4, 



C 2,J. - 75 (-fo,13 - ^0,14) - 75-^13,26 - 75-Pl4,26 (B.4) 
— \ (-Pll,21 — P.1,22) — 75 (&-Pll,23 + &P\\,14) 
~ \ (-Pl2,21 + -Pl2,22) + 75 (tt-Pl2,23 ~~ bPi2 t 24) 
+ | (-^21,31 + ^21,32) + \ (-^22,31 — -P22,32) 
+ 75 (a-P23,31 — W-23,32) — 75 (^24,31 + «-?24,32) 
+ 75 (^25,33 + ^25,34) ~ 75^33,4 + 75^34,4 
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where 



P a ,p = \E a ){E fi \ (B.5) 

and 



a = ai + (i2 
b = di — a 2 . 



(B.6) 



To obtain the annihilation operators of the I— th dimer, the index / should 
be added (c 1(2))tr — > c J)1(2)i<T , P Q;/3 -> P^) in (B.1)-(B.4). 
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Appendix C 



To present the Hamiltonian H in a compact form we first introduce the 
following operators (fj,, v — 1, 2) : 



S I,ti;J,v — a J,u,l ai ^h 
a 

S I,V,J," = \ ( n I,V,J,v,1 ~ n I,fJ.;J,v;i) ; 



a 

§J,H-,J,v = \ {lkl,ii.;J,u;l ~ Hl,,j,;J,u;l) 



=I,fj,;J,u =J,u,i 



Ti — , Ti 

— k [n T 7 „ — n T T , ) 

2 \=^,M;J,f;T =i,ij>;J,v;IJ 



(C.l) 
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The operators s z /^. Jl/ , s z /^. Jv , s^.j^, are not strictly the spin operators, 
they, however, show some similarities to the true spin operators (as e.g. (11)- 
(13)) and therefore they can be called "hopping spin operators". Moreover, 
the following abreviations have to be used: 

Q (I, fi; J, v) = ~3 Itli ■ (~£ I , mv + , (C2) 

R (J, fi; J, v) =Si tll ■ Sj, v ni,nJ,v, (C3) 
(/, n J, v) = (S'i^ - Sf^S'jJ « W (C4) 

R~> + (/, M ; J, „) = {S^S+ U - S+^SjJ *W (C5) 

The Hamiltonian if (cf (24)) with the use of (C.l) - (C.5) takes on the 
form 

J 1,(7 

-20 E [<9 l; / + 1, 1) + <5 (/ + 1, 2; /, 2)] 

+2t5£[Q(/,2;/ + l,l) + Q(/ + l,l;/,2)] 
£ [i? (/, 1; / + 1, 2) + A*." (/, 1; / + 1, 2) 

+R g >+ (I, 1; / + 1, 2) + i?->+ (/, 1; / + 1, 2) 
+R (I + 1, 2; /, 1) + R z >~ (/ + 1, 2; /, 1) 
+R*> + (/ + 1, 2; /, 1) + iT>+ (/ + 1, 2; 7, 1)] 
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-t<? E [R (1, 2; / + 1, 1) + R z >~ (/, 2; J + 1, 1) 

+ J R 2 <+ (/, 2; / + 1, 1) + A" + (/, 2; / + 1, 1) 
+R (I + 1, 1; /, 2) + i2* " (/ + 1, 1; /, 2) 
+R Z > + (/ + 1, 1; /, 2) + iT>+ (7 + 1, 1; I, 2)] 

+t(35 EE [i? (J, /i; / + i, n) + R z '~ (i, n) 

I 11=1 

+R Z > + (I,fi;I+l,fi) + R~' + (I,n;I+l, fi) 
+R (I + 1, fi; I, n) + R*~ (1 + 1,(1-1, n) 
+R Z ' + (/ + 1, fi; I, fi) + R-> + (I + l,n;I, fi)} 



z 

-* E J2 at «, 



(C.6) 



-2t/3 E EQ(/,^J,7i) + 2t5E EQ(/,//;J,/x) 
-t ((3 2 + 5 2 ) E E ^; 4 /4 + a*; 4 /*) 

+R z ' + {I,fi;J,fi) + R-'+{I,fi;J,fi)\ 

+2tps E E [R (i, v\ >h + R z '~ (i, v\ I -P) 

i^j n=i 

(I, fi; J, ~p) + R- + (I, fi; J, JI)] 

where h\ is given by (14). It should also be noted that the operators a T 
(o^ )) entering in (C.l) and (C.2), are in fact, dependent on (3 and 5 (see 

(29)). The decomposition of the Hamiltonian H (C.6) according to the terms 

23 



proportional to (3, S, (3 2 ,5 2 and (35 has thus only a formal character in order 
to keep the presentation of the Hamiltonian H in a compact form. 
The formula (C.6) is the complete expression. 
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